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ASYMPTOTIC SHAPE FOR THE CHEMICAL DISTANCE AND 
FIRST-PASSAGE PERCOLATION IN RANDOM ENVIRONMENT 

Abstract. The aim of this paper is to generalize the well-known asymptotic 
shape result for first-passage percolation on Z d to first-passage percolation on 
a random environment given by the infinite cluster of a supercritical Bernoulli 
percolation model. We prove the convergence of the renormalized set of wet 
points to a deterministic shape that does not depend on the random environ- 
ment. As a special case of the previous result, we obtain an asymptotic shape 
theorem for the chemical distance in supercritical Bernoulli percolation. We 
also prove a flat edge result. Some various examples are also given. 

First-passage percolation was introduced in 1965 by Hammersley and Welsh HW65 
as a model for the spread of a fluid in a porous medium. To each edge of the Z d 
lattice is attached a nonnegative random variable which corresponds to the travel 
time needed by the fluid to cross the edge. 

When the passage times are independent identically distributed variables, Cox 
and Durrett showed that, under some moment conditions, the renormalized set of 
wet points at time t almost surely converges to a deterministic asymptotic shape. 
Deriennic |Kes86| . and next Boivin |Boi90j . progressively extended the result to the 
stationary ergodic case. 

In this paper, we want to study the analogous problem of spread of a fluid in 
a more complex medium. On the one hand, an edge can either be open or closed 
according to the local properties of the medium - e.g. according to the absence 
or the presence of non-porous particles. In other words, the Z d lattice is replaced 
by a random environment given by the infinite cluster of a super-critical Bernoulli 
percolation model. On the other hand, as in the classical model, a random passage 
time is attached to each open edge. This random time corresponds to the local 
porosity of the medium - e.g. the density of the porous phase. Thus, our model 
can be seen as a combination between classical Bernoulli percolation and stationary 
first-passage percolation. 

Our goal is to prove in this context the convergence of the renormalized set of wet 
points to a deterministic shape that does not depend on the random environment. 
As a special case of the previous result, we obtain an asymptotic shape theorem for 
the chemical distance in supercritical Bernoulli percolation. 

In the first section, we give some notations adapted to this problem, precise the 
assumptions on the passage time and introduce the background of ergodic theory 
needed for the main proof. In section 2, we give some estimations on the chemical 
distance and the travel times. In section 3, we prove the existence of an asymptotic 
speed in a given direction and study its properties ( e.g. continuity, homogeneity 
and sub-additivity. . . ) Section 4 is devoted to the question of the positivity of 
these speeds, and in section 5 we prove the asymptotic shape result. In section 6, 
we study the existence of a flat edge in the asymptotic shape. In the last section, 
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we develop some examples, such as exponential times, correlated chi-square passage 
times, and a model for road networks. 



1. NOTATIONS, DEFINITIONS AND PRELIMINARY LEMMAS 

The percolation structure. Let us first construct a Bernoulli percolation struc- 
ture on Z d . Consider the graph whose vertices are the points of Z d , and put a 
non-oriented edge between each pair (x, y) of points in 7L d such that the Euclidean 
distance between x and y is equal to 1. Two such points are called neighbours. This 
set of edges is denoted by E d . 

Set n E = {0, 1} E . In the whole paper, p is supposed to satisfy to 

P G {Pc 1], 

where p c = p c {d) is the critical threshold for bond percolation on Z d . We denote 
by F p the product probability (p6i + (1 - p)8 )® Ed ■ 

A point oj in £Ie is a random environment for the first passage percolation. 
An edge e G K d is said to be open in the environment u> if w(e) = 1, and closed 
otherwise. The states of the different edges are thus independent under P p . 

A path is a sequence 7 = (xi, ei, X2, ei, ■ ■ ■ , x n , e„, x n +{) such that xi and a^+i 
are neighbours and e% is the edge between Xi and Xj+i. We will also sometimes 
describe 7 only by the vertices it visits 7 = (xi,x 2 , ■ ■ ■ ,x n ,x n+ \) or by its edges 
7 = (ei, 62, . . . , e„). The number n of edges in 7 is called the length of 7 and is 
denoted by |7|. Moreover, we will only consider simple paths for which the visited 
vertices are all distinct. A path is said to be open in the environment u) if all its 
edges are open in ui. 

The clusters of a environment u) are the connected components of the graph 
induced on Z d by the open edges in u). For x in Z d , we denote by C(x) the cluster 
containing x. In other words, C{x) is the set of points in Z d that are linked to x by 
an open path. We note x <-> y to signify that x and y belong to the same cluster. 
For p > p c , there exists almost surely one and only one infinite cluster. We denote 
by Coo the random set: = {H Z d ; |C(fc)| = +00}. Note that Coo is almost 
surely connected. 

Since p > p c , we have P p (0 G Coo) > 0. In our first-passage percolation context, 
it is natural to work on the infinite cluster. Thus we introduce the conditioned 
probability measure ¥ p on Qe by 

- P p Q4n{0eCoo}) 
p( } P P (0 e Coo) " 



The passage times. We give now to each edge of the graph Z a passage time, 
which represents the time needed by a fluid to cross the edge. Consider thus 
il s = R E , and a probability measure S v on fig . We use the abusive notation S v to 

help the reader to keep in mind the case of a product measure S v = v m , although 
we will consider the weaker assumption that £„ is a stationary ergodic probability 
measure - see below for a precise definition. 



The probability space. Our probability space will then be £1 = He x Og. A 
point in fi will be denoted (u),r]), with w corresponding to the environment and r\ 
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assigning to each edge a passage time. On O, we consider the probability P = ¥ p (E)S^ 
and its conditioned version: 

¥ = ¥ P ®S V . 

The two distances. We first introduce the chemical distance D(x,y)(ui) between 
x and y in Z d , depending only on the Bernoulli percolation structure u>: 

D {x,y){w) = inf |7|, 
7 

where the infimum is taken on the set of paths whose extremities are x and y and 
that are open in the environment to. It is of course only defined when x and y are in 
the same percolation cluster. Otherwise, we set by convention D(x, y) = +00. The 
distance D(x, y) is thus, when it is finite, the minimal number of open edges needed 
to link x and y in the environment to. We denote by j(x, y) an open path between x 
and y with this minimal number of edges. When necessary, we can define uniquely 
j(x, y) by choosing on order on the set of edges E d and taking the minimal such 
path for the lexical order on the edges of the path. 

For (lu, 77) G and (x, y) G Z d x Z d , we define the travel time d(x, y)(uj, rf) to be 

infVV, 

where the infimum is taken on the set of paths whose extremities are x and y 
and that are open in the environment lu. Of course d(x,y) = +00 if and only 
if D(k 1 l) = +00. It follows that the event {d(k, I) < +00} only depends on the 
Bernoulli percolation structure to. For t > 0, we note 

B t = {k G Z d ;d(0,k) < t}. 

It is the random set of points that can be reached from the origin in a time smaller 
than t, and it is a non-decreasing random set in t. 

Usual definitions and results in ergodic theory. We recall here some classical 
definitions, notations and results in ergodic theory which can be found in any course 
- see for instance the book written by Brown |Bro76| . 

Consider a probability space (O, J 7 , Q). A map 9 : fl — ► is said to be a measure- 
preserving transformation if for each A G T, Q(0^ 1 (A)) = Q(A). The quadruple 
(fi, J 7 , Q, 9) is then called a dynamical system. 

For any measurable set A, we can define the entrance time in A: 

n A (x) = inf{n > l;6 n x G A}. 

The well-known recurrence theorem of Poincare ensures that for Q almost every 
to G A, {n > 1; 9 n x e A} in infinite. 

If every subset of fl which is invariant by 9 has probability or 1, we say that 
9 is ergodic for Q. When 9 is ergodic, the mean entrance time is given by Kac's 
formula: 

riA dQ = 1. 



A 

For a set A G T such that Q(A) > 0, the entrance time in A induces a transfor- 
mation 9a '■ ^ — ► ^ defined by 9a{u) — 9 71a ^(uj). Of course, 9a is in general only 
defined for almost every lo £ A - this is a consequence of the recurrence theorem 
- but its definition may be extended to A. If Qa is defined by Qa{B) = Q< q^^ , 
the following results hold: 
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• If 9 is a Q-preserving transformation, then 9a is a Q^-preserving transfor- 
mation. 

• If 9 is ergodic for Q, then 9a is ergodic for Qa- 

Note that Kac's formula can here be rewritten in the following form: 

n A dQ A 



n 



Q(Ay 

The translation operators. 

Definition 1.1. For any set X , and u G Z rf , we define the translation operator 9 U 
on X E by the relation 

Ve G E d (9 u uj) e = u u . e , 

where u.e denotes the natural action of Z d on K d : if e — {a, b}, then u.e = {a + 
u,b + u}. 

In the whole paper, the measure S v which describes the passage time is supposed 
to be an ergodic measure preserving the Z d actions (9 u ) ul z%d . It means that for each 
u G Z d , (fis, B(Qs)> S U ,9 U ) is an ergodic dynamical system. This covers of course 
the case of a product measure = . 

Since is invariant under the action of 1 d , the law of the passage time of a given 
edge only depends from its direction. Let us denote by (ei, . . . , e^) the canonical 
basis of Z d . In the whole paper, we will suppose that 



(H int ) m= sup / V(0 , ej) dS u (r)) = sup Ve dS u (r]) < + x 

l<i<d J eSE d ' 

When needed, we will suppose that satisfies to 

(H a ) 3A, B > such that VA C E d S„ I ?/ G n s ; Vt)i> B|A| I < 



\ eGA 



A 

W 

for an appropriate value of a. Assumption (H a ) is in fact a combination between 
assumptions on the moments and assumptions on the dependence between the 
passage times of distinct edges. 

For instance, if S v is the product measure , assumption (H a ) follows from 
the Marcinkiewicz-Zygmund inequality as soon as the passage time of an edge has 
a moment of order 2a - see e.g. Theorem 3.7.8 in |Sto74| . Note that obviously, 
(H a ) is always fulfilled when the passage times are bounded. 

Definition 1.2. For each u G Z d \{0}, let T m (uj) — inf{n > 1; nu G Coo(oj)} and 
define the random translation operator on Q = Qe X &>s 

Note that T u only depends on the environment u>, and not on the passage times 
rj, whereas the operator Q u acts on the whole configuration (w, rj). The new config- 
uration Q u (uj, r]) is the initial configuration viewed from the first point of the form 
nu to be in the infinite cluster. 

Lemma 1.3. U is a V-preserving transformation. Moreover, O u is ergodic for P 
and 
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Proof. Let u G "L d . The dynamical system (f2, £>(£!), P, 6* u ) is the direct product of 
the mixing system (He, B(He), P p , 0„) by the ergodic system (O5, S(fis), S^, 
Since the product of a weakly mixing system and an ergodic system is ergodic - 
see for instance proposition 1.6 of Brown |Bro76| - , it follows that 9 U is ergodic 
for P. It is easy to see that Q M is the transformation induced by 9 U on the set 
A = {0 G Coo} x H s . Since ¥(A) = P p (0 G C 00 )5 I/ (fi s ) = P p (0 G <?«,) > 0, the 
announced result follows from the classical results of ergodic theory. □ 



Definition 1.4. Let us define for n > 1 

n-1 

T n>u H = ^T u (e^). 

fc=0 



Remark that T n>u represents the position of the n th intersection of the infinite 
cluster with the half line Z* + u. 

T nu 1 — 

Lemma 1.5. — 



5 P (0 G Coo) 



Proof. Since Q u is ergodic for P, this immediately follows from the pointwise ergodic 
theorem. □ 



2. ESTIMATIONS ON THE CHEMICAL DISTANCE AND TRAVEL TIMES 

We first recall the fundamental result of Antal and Pisztora (Theorem 1.1 in 
AP96 ): If p > p c (d), there exist a constant p = p(p, d) 6 [1, +00) and two strictly 
positive constants A and B such that: 

(1) Vy eZ d , F p (0^y, D(0, y) > pWy^) < Aexp (-I%[[i) . 

Note that if y = (y x , . .. ,y<i) G M d , \\y\\i is defined as follows: \\y\\i = Y,i=i \Vi\- 

Lemma 2.1. If (H a ) holds with two positive constants A and B, and if p is the 
constant given by £7J), then there exists a positive constant A' such that 

V.t eZ d , V(0<->x;d(p,x)>Bp\\x\\i) < j^—. 

Fill 

Proof. Let x G Z d . By definition, d(0,x) < J2ee-y(o x) 7 ?( e )- Then 

P(0 ^ x;d(0,x) > Bp||a;||i < P O^x; ^ 77(e) > Bp||ar||i 
We can now cut this probability in two pieces: 
O^x; ^ 77(e) > SpHxlU 

ee7(0,x) 

< P ( Q^x;D(0,x) <p||x||i; ^ rj(e) > £p||a:||i ) +P(0 <-► a;;D(0,ar) > p 

e67(0,a:) 
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Note that we already know that the second term decreases exponentially fast thanks 
to 0); for the first one, we have: 



0~x]D{0,x)<p\\x\\i; v(e)>B P \\x\\ 1 

7.|7|<P|MU V e 67 / 

P P (0 - «; 7(0,1) = y)S u ($>(e) > Bp\\x\\A 

y|<p||a;||i \eS7 / 

< Yl p p(° ~ =7)&,(l>(e) >fl| 7 |) 

7>l7l<PlW|i \ee7 / 



7>|7|<P|| 



< ^ P p (0 «-> x; 7(0, 1) = 7) F ^ with Hypothesis (JJ a ) 



7>l7l<PIM| 



AP„(0^a;) 1 . . . 
< T na as | 7 (0,x)| > ||x|j 



□ 



Lemma 2.2. There exists two positive constants A' , B' , depending only on the 
dimension d, such that for every r > 0, for every x £ Z d such that \\x\\i < r, 

P p (0 <-> x; D(0, as) > (3r) d ) < A' exp (—B'r) . 
Proof. Take p, A, B as given by (JTJ. Assume that r > max < 2, — j- I >. Take 



d-l 

.x g Z d such that ||a;||i < r. 

We denote by -B(r) the cube [— 5r/4, 5r/4] d n Z d . As the segment [a, 6] contains 
at most b - a + 1 integer points, |B(r)| < (^ + and, as r > 2, |B(r)| < (3r) rf . 
Denote by dB(r) the external boundary of B(r), defined as the set of points in 
Z d n B(r) c that admit a neighbour in B(r). Note that x G -B(r) and that if 
yedB(r), then Hi-yHx > |. 

Assume that D(0,x) > (3r) d , then D(0,x) > \B(r)\, and then 7(0,2:) must visit 
at least one point in dB(r): 

P(0 <-> x, D(0,x) > (3r) d ) 

< P ( 3y G &B(r) such that y, y <-> x, 13(0, y) > ^y- or £>(x, y) > ^y- 

< J] P /'o«j / , J D(0,y)>^ N ) +P ^ y>jD ( :r! y)>^rLy 

yedB(r) ^ ' ^ ' 
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(3r) d 

But — - — > bpr > pmax{||a;||i, ||x — y||i}, and thus 

P(0 <-> x, D(0,ar) > (3r) d ) 
< Y, W(0~y,D(0,y)>p\\x\\ 1 )+W(x~y,D(x,y)>p\\x-y\\ 1 ) 

yedB(r) 

But min{||x||i, \\x — y\\{\ > § and using (IJ we get: 

P(0 <-> x, L>(0,x) > (3r) d ) < 2 ^ Aexp (-- ^ J 

< 2(2d)(5r/2 + l) d - 1 J 4exp 
This ends the proof of the lemma. □ 

Lemma 2.3. If (H a ) holds, then there exist positive constants a and A' , depending 
only of the dimension d of the grid, such that for every r > 0, for every x G Z d 
satisfying \\x\\i < r, 

A' 

P(0 <-» x,d(0,x) > ar d ) < —, 

with v = a — d + 1 . 

Proof. Let b > be a constant that will be fixed later. Denote by B r the cube 
[-(1 + b)r, (1 + b)r] d n and by ext .B r the set of points in Z d \B r that admit a 
neighbour inside B r . We have 

|-B r | < (2(l + 6)r + l) rf . 

Denote by n e (B r ) the number of edges thqt have their two extremities in B r . Then 

d(2(l + 6)r - 2) d < n e {B r ) < d(2{l + b)r) d . 

Now, let A and B be the two constants given by Hypothesis (H a ), and choose a > 
such that 

Vr > 0, a d > Bd(2{l + b)r) d . 
Take x G Z d such that |jx||i < r. Note that x G £?„. Now, 
P(0 <-> x, d(0,x) > ar d ) 

= P ( <-» x, d(0, x) > ar d and ^ 7j e < Bd(2(l + b)r) d J 

+P ( <-> x, d(0, x) > ar d and r/ e > Sd(2(l + fe)r) d j 
< P I By G d ext B n such that <-> y, y <-> x, d(0, y) > — or d(x, y) > — 



+P ( - Bn ^ B r)j 

\e£B r J 
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By Hypothesis (H a ), the last term is smaller than — ; — — . Note moreover that 
v ; n e (B r ) a 

for r large enough: 

Bp\\y\\i < Bp{(\ + b)r + 1) < ar d /2 and \\y\\i > br 

Bp\\x-y\\i<Bp((l + b)r + l)<ar d /2 and \\y-x\\i>br. 

Thus 

3y 6 d ext B r such that <-> y, y <-> x, d(0, y) > — or d(x, y) > — 

< Y p(o<-»y,d(0,tf) > +p(x^y,d(x iy ) > 
yed ext B r \ * J \ * 

< Y, y(0^y,d(0 1 y)>Bp\\y\\ 1 )+¥(x^y,d(x,y)>B P \\x-y\\ 1 ) 

- A( >T1& + II 1 11° ) with Lemma EU 

< 2d(2(l + 6)r + l) d - 1 ^. 

□ 

Lemma 2.4. There exist a function ip : 7L d — > K + satisfying \im^ u ^ 1 ^ +OQ ip(u) = 0, 
and suc/i i/iai 

V W £Z d , E gp OP(0,7») < Pp( ^ 1 oo) (P + ^))- 
Proof. Remember the following identity: 

E f ^Tir^= / °°f P (D(0,T u u)>h\\u\\)dh, 

" \\ u \\ JQ 

and cut E p (-D(0, T u u)) in three terms, according to the values of T u . Let us 
estimate first the contribution for T u small: 



Si(h,u) = p P fi>(o,r T ,u)>fcHi,r u < 



3||«[[i 



1 


P p (0 <-i 


• 00 ) 


1 




P p (0 <-> 


• 00 ) 


1 




Pp(0 4-) 


• 00 ) 



»p(i?(o,r t ,u)>/i||u||i,r 1 ,< 



3||u||i 

Y P P (0^ ku,D(0,ku) > ft||u||i,T u = fc) 



- F_m ^ ^ £ Pp(0~fcu,I?(0,AtO>fc||u||i). 

■ (feii"iii) 1/d 
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(7i||tt|| 1 ) 1 / d 

But as k < n — n , we have /illulli > (3fc||u||i) d and thus, by the previous 

3||u||i 

lemma, 



SliKu) ~ P (0 L oo) £ F p (0^ku,D(0,ku)>(3k\\u\\ 1 ) d ) 
S P p (0-oo) 3||«||i 6P v 3 J' 



Now, J R+ S\(h,u)dh is finite, and there exists a constant Ci, independent from u, 
such that 



f p {Q^oo)J 3||u||i V 3 / 

< i a' r°° , IA ( „,xW 







< 



P p (0 <-> oo) 3||u|| 2 
Ci 



Let us now estimate the contribution for T u medium: 



■*2i/'.«) = IF' ( D(0,T u u) > hWuh, {k \^f /d < r M < - 

3||u||i p 



> P I £>(0,T u u) > feUuHx, {h \ U ^ 1/d <T U <- 



P p (0^oo) "V 3|M|! a ~p 

= p f0 ^ x E Pp(i>(o,r u «)>/i|H| 1 ,r u = fc) 

31R1R" ^ - p 

^ iTTTT V E P p (0^fcu,D(0,A5tt)>/t||tt||i) 

P o (0 <-> oo) ' 

^ iTTTT T E Pp(0~to,I>(0,*u)>P*Hi) 

P B (0 <-> oo) • <: — ' 

and with (JIJ: 

^ iTTTT V E Aexp(-Bfc|| u |j 1 ) 

P D (0 <-> oo) ^ 

~ P p (0 ^ oo) p ^ V 3 
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Thus, there exists a constant C%, independent from u, such that 

S 2 (h,u)dh < — [°° -Aexp(-B ^ U ^ 1/d ) ,//, 

< ^ _A_ [°° x Vd exp f -S^) dx 



< 

Finally, for T u large: 



c 2 



and 



■S'. : (/'.'0 = IP,, (D(O i r,u)>%|i l T 4 > ~ 
< P P [T U >- 



J S 3 (h,u)dh < J ¥ p (r u >^jdh 

< pE p T u with Lemma 1 1.31 



? p (0 <^ oo)' 

Putting these three pieces together, we get the announced result. 

Corollary 2.5. For every u <E Z d \{0}, 

E ¥ (d(0,T u u)) < mE Fp (D(0,T u u)) < oo. 

Proof. 

oo 

E F (d(0,l») = ^E F (d(0,/c)l {Tu=fc} ) 

fe=0 

OO / 

fc=0 y eS7(0,feu) 

fc=0 \ \e£7(0,fcu) 



□ 



< J2 E p p {MT u =k}mD(0,ku)) 

fc=0 

< mE ?p (D(0,T u u)). 

The result of the previous lemma gives then the integrability of d(0, T u u). □ 

3. Existence of directional speed constants and related properties 

The next Lemma is the analogous of the first step in the classical study of first- 
passage percolation, that is the proof of the existence of an asymptotic speed in a 
given direction. 
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Unsurprisingly, the fundamental tool is here again Birkhof 's subadditive ergodic 
theorem. The main difference is that it is applied to an induced dynamical system, 
and not directly. It is the reason why checking the integrability flemma l2.5fl was 
more intricate than in the classical case. 

Lemma 3.1. For n > and u 6 Z rf \{0}, we define the travel time f n . u between 
the origin and the n-th intersection of the infinite cluster with Jf+u: 

/n,«(w,?7) = d(0, (T ntU (u>)u)(uj,r)). 

Then there exists a constant f u > such that 

^ -» /„ P a.s. 
n 

The convergence also holds in L 1 (P). Moreover, f u < E pd(0, T u u). 

Proof. For the convenience of the reader, we now choose u € Z d \{0}, and we write 
/„ (resp. /, T n , T) instead of (resp. f u , T n , u , T u ). 

Note that the travel time is a sub-additive function, we thus have: 

fn+k = d(0, T n+k u) 

< d(0, T n u) + d(T n u, T n+k u) 
= d(0,T n u) + d(0,T k u) oT n 

= fn + fk°T n . 

Moreover, the /„ are non-negative and, thanks to lemma fOl P is invariant under 
<d u . Then, as f\ is intcgrable thanks to lemma Birkhof 's sub-additive ergodic 
theorem says there exists a O u -invariant function / such that 

— / P a.s. and in L^P). 
n 

Since <d u is ergodic for P, every u -invariant function is P almost surely constant, 
so / is constant. Moreover, the deterministic sub-additive ergodic theorem ensures 
that 

E -f 

Ei/= inf < E ™/i = E ¥ d(0, T u ). 

n>0 n 

□ 

We can now define a function /i which is exactly the analogous of the function 
that appears in the classical case. 

Theorem 3.2. Suppose that (H a ) holds for some a > 1. For each u S Z d \{0}, we 
set n(u) = P p (0 <-> oo)/ u , where f u is given by the previous lemma, and /i(0) = 0. 
Then 

,. d(0,T nu u) - 
hm ■ = filu) P a.s. 

"- > °° T n , u 

Moreover, this function /i : 7L d — > [0, oo) has the following properties: 

(1) /i is symmetric: for any vector u £ Z d ; fi(—u) — n(u). 

If, moreover, Si, is a product measure, then \x is invariant under permuta- 
tions or reflections of the coordinates in Z d . 

(2) /i is homogeneous, in the following sense: for any integer k G Z and any 
vector u £ 1 d : 

/j,(ku) — \k\fi(u). 



12 



CHEMICAL DISTANCE AND FIRST-PASSAGE PERCOLATION 



(3) p, is bounded, in the following sense: for any vector u 6 Z d : 

< pTi||w||i and p(u) < /i*||w||i, 

with /j* = sup 1 < 4 < d /j(e 4 ). 

(4) /i is sub-additive, in the following sense: for any vectors u and v in1 d , 

p(u + v) < p,(u) + n(v). 

(5) p is continuous, in the following sense: for any vectors u and v in1 d , 

\p(u) - p{v)\ < - v\\i. 
We delay the proof of this theorem to make some remarks and comments. 
Remarks. 

• We exactly recover the properties that the results of classical first-passage 
percolation let expect. 

• Let us now discuss the meaning of the inequality p(u) < /9m||u||i. In 
classical first-passage percolation, we would have the inequality p(u) < 
m||u||i. Roughly speaking, we can say that the additional factor p = p(p, d) 
corresponds to the bound given by corollary 1.3 of Antal and Pisztora^for 
the asymptotic ratio between the chemical distance and the l\ distance on 
Z d : 

(2) lim" ^nr^tf"*} < P(P, d) P P a.s. 

• In the special case S v = (<5i)® E , the travel times coincide with the chemical 
distances. As a direct corollary of the previous theorem, we obtain the 
existence of asymptotic directed speeds for the chemical distance: For each 
y G Z d \{0}, there exists a constant p(y) such that 

D(0,ny) 
lim v yj = fi(y) F p a.s. 

0<-my 

Since p(y) < p(p, d)\\y\\i, this precises the equation © of Antal and Pisz- 
tora. 



Proof. Step 1. The convergence result is clear since T n , u /n — > 1/P P (0 G Coo) and 
fn,u/n — * fu > 0. Consequently, we obtain the identity fi(u) = P p (0 <-» oo)f u . 

Step 2. It is easy to see that d ^°'^ u ^> weakly converges to 

P p (0 «-» oo)S fu + (1 - P p (0 «-► oo))6 +oc . 

By translation invariance d(0, nu) = d(nu,0) has the same law that d(Q, — nu). It 
follows that f u = f- u , and then that p(u) = p(—u). 

When S u is a product measure, Property (1) is obvious since the whole model is 
invariant under permutations or reflections of the coordinates in Z d . 
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Step 3. Let us prove now that for any integer k G Z\{0} and any non null vector 
u e Z d , fj,(ku) = \k\p(u). Note that 



fn.ku I , I fn.ku 

- \k\ 



n n\k\ ' 

but (f n ,ku/n\k\)n>i is a subsequence of {fn,eu/n) n >i with e = k/\k\. Thus 
lim '^" ,fc " = |fc| lim '^" ,£ " = |fc|/i(eu) = |fc|/i(u), 

n— »+oc n n^+oo n 

by the previous step. 

Step 4. By the previous step p(u) = -^p(nu). Then 

u(u) = P(0^oo) — 
n 

E ¥ (d{0,T nu nu)- 

< — F(0 <-» oo) 

n 

mE w (D(0,T nu nu) -, n 

< - F(0 <-> oo) 

n 

< Mp + ^Wnuh by lemma m 

< m(p + ip(nu)\\u\\i 

Since lim n _> +00 ip(nu) = 0, it follows that < mp||ti||i. This proves the first 

part of (3). 

Step 5. Let us now prove the sub-additive property of fi. Let u and v be two fixed 
non null vectors in Z d . The main difficulty here, compared to the classical proof in 
standard first-passage percolation (see Kesten) is that two given points in 1 d are 
not linked by an open path in every configuration u>. 

Lemma 3.3. Let x and y be two linearly independent vectors in "L d . Then 
1 - 

— l{fca:^oo}l{fcanoo} Pp(0 <-> oo) 2 almost surely under ¥ p . 

k=l 

Proof. Let a = | min(||x||i, ||y||i, ||a;-J/||i) and define Z k = TL {ak < c(kx)} l {ak < c{ky)} . 
P p (Z k ^ l {fccwco} l {fe „ MOo} ) < P p (afc < \C(kx)\ < +oo) +P p (afc < \C(ky)\ < +oo) 

< 2¥ p {ak < \C(0)\ < +oo). 
Since p > p c , C(0)l{|c(o)|<+oo} is integrable, and we get 

^2Pp(Z k ^ l{fcr^oo}l{/cywoo}) < +00- 
k>l 

Hence, Borel-Cantelli lemma ensures that, P p almost surely, Z k coincides with 
lL{kx^oc}TL{ky*->oo} for every k large enough. It follows that ^J2l=i l{fei^oo}l{fcy^oo} 
has the same behavior that ^ X)fc=i But (Z k ) k >i is a sequence of independent 
uniformly bounded random variables. Then 



1 n 



n 
fe=i 
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Since a = § min(||x||i, ||y||i, ||x - y\\i), 

E Pp Z k = E Pp l{ak<\C(kx)\}K F p l{ak<\C(ky)\} = V p {ak < C(0)) 2 -> P p (0 <-> C50) 2 , 

and the results follows from Cesaro's theorem. □ 



Lemma 3.4. Suppose that there exists a > 1 such that (H a ) is fulfilled with two 
positive constants A and B and let p be the constant given by 0J). Let x and y be 
two linearly independent vectors in 2i d , and for k E denote by A k the following 
event: 

A k = {d(0,kx) < Bpfc||a;||i, d(0, k(x+y)) < Bpk\\x+y\\i, d{kx,k{x+y)) < Bpk\\y\\i}. 
Then, almost surely under P, if 

1 ^ d(0,kx). 



-J2 

n ± — ' 



n * — ' k 

fc=l 



-1a,1i 



A k - 1L {0^oo}l{fcx^oo}- IL {fc(i:+j/)^oo}j 



I, 



then (u n )n converges to /i(x)l{owoo}Pp(0 oo) 2 and moreover, (E p p (u n ))n con- 
verges to /i(a;)P p (0 <-» oo) 3 . 

Proof. Note first that 

}^{k(x+y)^oo} # 1{0^oo}1{ kx^oo}^{k(x+y)^ce}*) 
= E Pp (l{o^oo}l{fea;noo}l{fc(x+ a )^oo}(l ~ IaJ) 



< 



\^-{Q^co}^{kx^oo}^{k(x+y)^cc} 



x (l{d(0,fex)>Bpfc||x||i} + l{d(0,fe(a+y))>Spfc||x+y||i} + ^{d(kx,k(x+y))>Bpk\\y\\i})) 

< P p (d(0, fere) > B/9ft[[x||i and «-> fex) 

+P P (d(0, fc(x + y)) > Spfc||a; + and <-> fc(x + y)) 
+P P (d(fcx, fc(x + j/)) > and fcx <-» fc(x + j/)) 

A / 1 1 1 



< 



*° llWIf \\x + y\\? WvWi 



with the previous lemma. 



Then, using Borel-Cantelli Lemma, we know that almost surely there exists K = 
K(ui) > such that for any k > K, 

lA fc l{Owoo}l{fex^oo}l{fe(a;+jf)^oo} == l{0^oo} l{fcx^oo} l{k(x+y)<->oo} ■ 

Choose e > 0. Then, as a consequence of the convergence result in Theorem l3.2l 
and enlarging K = K(oj) if necessary, one has, almost surely: 



Vfe > K(w), 



d{0, kx) 



k 
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Now, for every n: 



\u n - /i(x)l {0 ^ oo} Pp(0 <-> oo) 2 

K(uj) - 1 d(o,k x )_ 



1 f d(0,kx) 
n K 



n 

+ 1 E 



k=K(ui) 



d(0, kx) 



k 



l^fcl{0^oo}l{fcK^oo}l{fe(x+i/)^oo} 



1 n 

- E ^^{O^oojlffcr^oojlffefir+^oo} - l{0«oo} P p(0 <-> OO)' 



k=K(uj) 



But we can choose 7V(u;) large enough to ensure that if n > N(u>), then the first 
term in the sum is smaller than s. The second one is smaller than e thanks to 
the choice of For the last term, remember that the choice of K(w) ensures 

that for every k > K(u>), we can forget the indicator function lA fc - But then the 
previous lemma ensures the almost sure convergence to of the last term. This 
proves the desired almost-sure convergence. 

To prove the convergence in mean, note that u n < p and use the dominated 
convergence theorem. □ 



Let us now prove the sub-additive property of p. Note first that if x and y are 
linearly dependent, the inequality is in fact an equality thanks the the homogeneity 
property. Now, let x and y be two linearly independent vectors in Z d . For every 
k > 1, let Ah be the same event as the one defined in the previous lemma. Then, 
for every k > 1, 



d(0,k(x + y)) 

^ R A k -"-{O^oo} -"-{fcir^oo} R {k(x+y)^oo} 

d(0,kx) 

7 *-A h JL{0^oo} ^{kx^oo} *-{k(x+y)^oo} 



< 



d{kx,k(x + y)) t 

T ^ JlA fe IL {0<-yoo} ,L {kx^oo} a -{k(x+y)<-^oo} 



Summing for 1 < k < n and taking the means gives: 



(3) 
(4) 
(5) 



E 



+ E 



'if dfl.fcQc + y)) , , , 



k=l 



(If d(0, for) 

< p k ^ A k^{0^co}^{kx^co}l{k(x+y)^oo} 



k=l 

1 f» d(kx, k(x + y)) 



E 

fe=i 



lA fc l{0^oo}l{fcir^oo} l{fe(a;+j/)^oo} 
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Thanks to the previous lemma, J2J and J3J converge respectively to /i(x + y)P p (0 
oo) 3 and /j,(x)¥ p (0 «-> oo) 3 . Moreover, 



1 ™ 

T? ' 



/ d{kx, k{x + y)) 



n f — ' V 



lA fc l{0^oo}l{fei;noo} l{fe(K+y)^oo} 



fc=l 



n k=l ^ 

where is the event 

B fe = MO, -ky) < P Bk\\y\\, d(0, -k{x+y)) < P Bk\\x+y\\, d(-ky, -k(x+y)) < P Bk\ 

This was obtained by translating along the vector — kx. Once again, the previous 
lemma enables us to conclude that (JSJ) converges to /z(y)P p (0 <-> oo) 3 . As P p (0 <-> 
oo) > 0, this leads to the desired inequality. 

Step 6. Property (5) is a direct consequence of the sub-additivity of \i. 

Step 7. Using properties (1), (2) and (4), we get: 

d 

fi(ui, . . . , u d ) < \v,i\ii(ei) 

i=l 

where (ei)i<i<d is the canonical basis of 7L d . Thus, 

The second part of property (3) follows immediately. □ 

This finally allows us to prove an analogue result for any rational direction. 
Indeed, let q £ Q d \{0} be fixed. Choose any couple (TV, u) £ Z* + x 7L d such that 
u = Nq £ 7L d and set fi(q) = p(u)/N. The constant /j,(q) is well defined thanks 
to the homogeneity property of fi. The function fi, extended in this manner to 
Q d \{0}, obviously keeps its properties: for any rational a and any vectors u and v 
in Q d , 

p(au) = |a|/i(u), and fi(u + v) < p(u) + p(v) and \p(u) — (i(v)\ < pm||u — v\\x. 
Thus we can extend /i to M. d by continuity, keeping these three properties. 



4. Sufficient conditions for the positivity of fi 

As in the classical cas, we study now the question of the positivity of /i, which is 
a crucial point in the quest of an asymptotic shape result. We begin with the case 
of a product measure. 

Theorem 4.1. Suppose here that S v — v E . Then 

• If pis (0) < p c , then fi is a norm on M. d . 

• If pv(Q) > p c , then \i — on R d . 

Proof. The only remaining point to consider is to determine whether p(u) = 
implies u = or not. Let us define, for each e £ K d , r]' e = n e uj e + (1 — ui e ). The 
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law of (?7e)eeE d 

under P is (pv + (1 - p)S 1 ) tg,E . We can define the new distance 
relatively to these new passage times: 



where the infinimum is taken over all the paths from x to y in Z d . This is in 
fact the classical first-passage percolation distance associated to the passage times 
v 1 = pv + (1 — p)5i. Denote by p! its associated renormalized limit. Note that by 
construction d'(x,y) < d(x,y) 7 and thus that p! < p. It is a now classical result 
(see Kesten |Kes86| for instance) that z/(0) < p c implies that the associated limit 
p' is a norm on M. d , and, since ^'(0) = pv{fS), if pv(0) < p c then p is a norm on M. d . 

On the other hand, suppose that pv(Q) > p c . Color in red the edges e such that 
Lu(e) — 1 and z/(e) = 0. As pv(Q) > p c , there is almost surely a unique infinite 
cluster of red edges, included in the infinite cluster for the Bernoulli(p) percolation 
structure. Let u G % d . On the event R = "the origin is in the red infinite cluster" 
(included in the event "the origin is in the Bernoulli(p) infinite cluster"), we can 
find an increasing sequence of (random) integers (k n ) such that k n .u is in the red 
infinite cluster for every n. But then clearly, d(0,k n .u) — 0. Thus, evaluating fi(u) 
on this subsequence gives fJ.(u) — on the event R. Since /x(m) is constant under 
P, this completes the proof. □ 

Corollary 4.2. The asymptotic speed fi associated to the chemical distance in 
Bernoulli percolation with parameter p > p c (d) is a norm on IR^. 

When S v is not a product measure, there is no general method to determinate 
whether v is a norm or not, even in the classical case p = 1. 

A natural idea is to use stochastic comparison: if S v i >- S v , then n(p,v') > 
Pl{jP,v). Consequently, if a is such that S v i(r) e S [a, +oo)) = 1 holds for each 
e G E d , then S v > is a norm for each p > p c , because d(., .) > aD(., .). 

We will now prove that p, is a norm if S v satisfies an appropriate large deviation 
inequality. 

Theorem 4.3. Let us suppose that there exists A > and a map f : (0, +oo) — > ]R+ 



Let K = lim exp(— /(e)). 

Then, the application p associated to S v is a norm as soon as Kop\(d) < 1, 
where X(d) is the d-dimensional connective constant 




Ve>0 VACE d & Uefi s ;$^>J*<e|A| < Aexp(-/(e)|A|). 



Proof. Let e > be small enough to ensure that (X(d) + e)exp(— /(e)) < 1 and 
consider u G Z d . 



18 



CHEMICAL DISTANCE AND FIRST-PASSAGE PERCOLATION 



P(d(O lU )<e||«||i) = ^P(7(0, U )-7;E^^ £ ll u Hi) 

7 e£A 

< ^P( 7 (0, U )= 7 ;^^<£|7l) 

7 eeA 

< H(X(d)+e) n exp(-nf(e)) 

™=ll«lli 

^ i TTU^T^ r -^T((A(d) + £ )cx P (-/(e))) 11 " 111 . 

1 - ((Ma) +e)exp(-/(ej) 

Then, it follows from Borcl-Cantclli lemma that d(0,nu) > n-|[ti|| j if n is large 
enough: we can conclude that fJ.(u) > e\\u\\i. As fi is homogeneous and continuous, 
this inequality is extended first to Q d , and next to M. d . Since ||.||i is a norm, so 
does p. 

□ 

Remarks. We will give in subsection !7 . 2l an example of application of this theorem. 

Let us now discuss the assumption of theorem 14.31 by considering the case S v = 
(Ber(l — q))® E . By Chernof's theorem, one can take 

/(e) = (1 - e) log — + slog- " 



q l-q 

thus Kq = q. Then, theorem 14 . 31 savs that p is a norm as soon as qp < whereas 
theorem 14. II teaches us that the optimal condition is qp < p c . 

Since ~ 53 ~ Pc(d) for large d, assumption of theorem 14.31 does not seem so 
bad. 

5. The asymptotic shape theorem 

From now on, we suppose that fj, is a norm. We want now study the convergence 
of the renormalized set of wet points at time t to the unit ball for fi. Because of 
the presence of holes in the infinite cluster, the Hausdorff distance seems naturaly 
adapted to study the convergence. 

Definition 5.1. For x € M. d and r > 0, we define 

B f ,(x,r) = {yeR d ;[i(x- y ) <r}. 

Definition 5.2. The Hausdorff distance between two non empty compact subsets 
A and B of M. d is defined by 

V(A, B) = inf{r > 0; A C B + B^O, r) and B cA + B^(0, r)}. 

Note that we use the Hausdorff distance associated to fx, but the equivalence of 
norms on R d ensures that the induced topology does not depend on this choice. 

Theorem 5.3 (asymptotic shape). Suppose that (H a ) holds for some a > d 2 +2d—l 
and that p, is a norm. Then, 

lim V ( — ,BJQ,1) ) = ¥ a.s. 
\ t J 
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Remark. It is important to have in mind that (H a ) is always satisfied if the 
passage times are bounded by an absolute constant. We have already seen that S v 
is always a norm when the passage times are bounded from below by a a positive 
constant. Putting these two results together, we get that S v always satisfies the 
assumptions of theorem 15.31 if it is an invariant ergodic measure for which we can 
found < a < b < +oo with S v ([a, b] E ) = 1. We thus obtain an asymptotic shape 
theorem for the chemical distance in supercritical Bernoulli percolation. Remember 
that on the event "the origin is in the infinite cluster" , we denote by B n the set of 
points in Z d whose chemical distance from is less or equal to n. 

Corollary 5.4. For every p > p c (d), there exists a deterministic convex compact 
set A with non-empty interior and invariant under permutations or reflections of 
the coordinates in M. d such that: 

► A for the Hausdorff topology P p a.s. 

n 

We will first give some preliminary lemmas in order to prove theorem 15. 31 
Lemma 5.5. Let z £ 7L d be a fixed non null vector. Then for every e > 0, 

¥(3N >0,VxE [N, +oo), 3k £ [(l-e)x, . . . , (l+e)x}r\Z + such that kz <-> oo) = 1. 

Proof. Almost surely under P, the sequence (T m , z ) m >i is an unbounded increasing 
sequence; thus there exists m £ Z+, with T m , z < x < T m +i, z . Now, 

X Tra^z 

which vanishes at infinity since Tm,z tends to w?r^ — It follows that x ~ Tm ' z < g 
as soon as x is large enough, thus we can take k = T m z . □ 

Lemma 5.6. Suppose that (H a ) holds with a > d 2 + 2d — I and let e > be small 
enough. Then 

P(3M > 0, Vm > M, Vx £ Z d ,Vy £ Z d , (|[ac||i = m, \\y\\i < em) => d(x,y) < Bpem) = 1. 
Proof. In view to use Borel-Cantelli lemma, denote by A m the set: 

A rn = {3x £ Z d , 3y £ Z d such that ||x||i = m, \\y\\i < em and d(x,y) > Bpem}. 

Using lemma 12731 and lemma I^Tl we get (note that in the following lines, Ki and Lj 
denote constants, depending only on the dimension d, whose values are not precised 



T, 



771+1,2 



1 = 
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to avoid intricate expressions). We set r = ( Bp ^ m ^ 



1/d 



V{A m ) < E E E F(x~y,d(x,y)>Bpem) 

x,||x||i=m 1=0 y : \\x—y\\i=l 
r 

^ E E E P(y -0,%,0) > Bpem) 

xJMI^m i=0 y,||y||i=( 
em 

+ E E E P(»^0,d(0,y)>Bp||i,||i) 

x,||x||i=m i=r j/,||t/||i=i 

< f y y — 

— Z-^l Z-^l Z-^l j.v 

x,||x||i=m Z=0 y,||j/||i=; 

em j 

+ E E E £ 

x,||x||i=m i=r |/,||i/||i=i 

Let us compute the first term 

x,||x||i=m 3/)l|y||i = ^ x,||x||i=m /— 

r d 

— ^ 2 / / 

x,||xj|i=m 



and the second term: 

em j em T 

y y y hi < y V;"-^ 

x,||x||i=ro l=r y,\\y\\-i=l x,||x||i=m l=r 

+oo 

E E^f 

ic,||ie|| i — m l=r 

E L 



< 



< 



m (a-d)/d 
x,\\x\\\—m 



< 



L 



3 



Remember that v = a — d+l. The choice we made on a ensures that both terms 
are the general terms of convergent series, which allows us to conclude thanks to 
Borel-Cantelli lemma. □ 

Lemma 5.7. Let e > 0. Then 

P(3M > 0, Vy e Z d , (||y||i > M and y |d(0, y) - < e||y||i)) = 1. 

Proof. Let us suppose this is false and choose e > such that the assertion fails. 
Then, with positive probability under P, there exists a random sequence (y n )n>o 
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of points in Z d such that 

\\Vn\\i -> oo 
Vn <->• 

\d(0,y n ) - fj,(y n )\ > e\\y n \\i 
By considering a subsequence if necessary, suppose that 

Vn _^ 
Il2/n||l 

Let us approximate z by a renormalized integer vector: consider e\ > (to be 
chosen small enough later), and choose z' 6 Z d be such that || y^ji — z\\i < E\. 

Let us find, for each y n , an integer point on the line M.z' close enough from y n : 
let h n be the integer part of ttt^-. We have 



\\y n -h n -z'h < 



Vn 



< WVnWl 



\\z'\\l 

Vn 



\\Vnh 



z' i 



/in 



Ik'lli 



IIjmIIi \\z'\\i 

Choose N > large enough to ensure that n > N => | 
With the choice we wade for z' we have n > N - 



z i 



J.V" 111 



.y« Hi 



z||l < £l, 

< 2^i, whence 



l z 'lli 



\\yn-K.z'\\ 1 <2e 1 \\y n \\ 1 + \\z'\\ 1 . 

Lemma 15.51 ensures now that, enlarging N if necessary, there exists for every 
n > N, a integer fc„ G [(1 — E\)h n , . . . , (1 + E\)h n \ such that kz' is in the infinite 
cluster. But: 

\\y n - k n .z'\\i < \\y n -h n .z'\\i + \h n -k n \\\z'\\i 

< 2e 1 [[y n [[ 1 +||/|| 1 +£ 1 /i n [[z'|| 1 

< 2e 1 || 2 /„|| 1 + ||z'|| 1 + £l ||y n || 1 

< 3ei||2/„||i + ||z'||i 

< 4ei||y„||i, 

if N is large enough. Enlarging once again N if necessary, we can use lemma 15.61 
to have: for every n > N: 

d{y n ,k n .z') < ABpsi\\y n \\i. 

Finally, for every n large enough, we have: 

\d(0, y n ) - n(y n )\ < \d(0, y n ) - d(0, k n .z')\ + \d(0, k n .z') - n{k n .z')\ + \fi{k n .z') - n(y 

d{0,k n .z') 



< ABpetWynW! + k n 



< 4Bpe l \\y n \\ l + (l + e 1 



+ (J,l\\k n .z' - y n \\l 



\\Vr, 



\* 1 



d(0,k n .z') 



fJ-(z') 



■4ei/ii||y n ||i 



But the convergence in the direction given by z' ensures that for every n large 
enough, 

d(0,k n .z') 



k n 



V(z') 



<ei. 



22 



CHEMICAL DISTANCE AND FIRST-PASSAGE PERCOLATION 



So by choosing e\ small enough, we can ensure that for every n large enough, 

\d(0,y n ) - n{y n )\ < e||2/™||i, 

which yields to a contradiction. □ 

Proof of theorem \b.ift Let e > 0. 
Step 1. If 

3T>0 Vt>T jCfl M (0,(l + e)) 

fails, then there exists an unbounded increasing sequence of times {t n )n>i and a 
sequence of vertices (y n )n>i with y n G B{t n ) and j 2 - ^ B^{Q, 1 + e), which gives 
d{0,y n ) < t n and /i(y„) > (1 + e)t n . By lemma l5"71 this happens with a null 
probability. 

Step 2. Let us show that with probability 1, the following property holds: 
(6) 3T>0 Vt>T S M (0,1-|) c ^+^(0,|). 

If © fails, then, with positive probability under P, there exists an unbounded 
increasing sequence of random times (t n )n>i and a sequence (v n ) n >i of points in 
R d such that 

f «„ € fl„(0, 1 - f ) 

\ v n i ^ + B„(0, f ) 
Let u be a limiting value for (u n ) n >i: clearly, there exists an unbounded increasing 
sequence of times t' n such that v £ -^- + 5^(0, |). Let w G Q d such that fx{v — w) < 

Then, for each n > 1, we have w ^ -p^ + £? M (0, |), whence it; G B M (0, 1 — ^e)- 
We can write w as w = z/Z, with z G Z d and Z G Z + . 

Lemma 15.51 ensures that, provided that n is large enough, there exists z n G 
Bjj,(t' n w, j[i(t' n w)) with 2 n «-> oo. The vertex z„ can not be in -jr 2 - because w £ 

^ + B^O, |). It follows that c£(0, z n ) > i„. But 

fi(z n ) < fi(t' n w) + -n(t' n w) < t' n n{w) + -t' n < t' n ( 1 - — + - J = (l - -J t' n . 

By lemma l5~7l this can (almost surely) not happen. Then, JfjJ) holds almost surely. 
Since 

B M (o,l-|) C + -B^ (o,|) =^B M (0,1) C^+S M (0,e), 
this completes the proof. □ 

6. Properties of the asymptotic shape 

As in the classical case, one has immediate properties for the asymptotic shape 
(see |Kes86| . chapter 6): 

Lemma 6.1 (Stochastic comparison). Let p c < p < p' and be two probability 
measures on R + such that the assumptions of the shape theorem are fulfilled for 
S v = v m and S v . = (*/)® E and such that v' is stochastically smaller than v, then 
if A and A' denote their respective associated asymptotic shapes, we have 

A c A'. 
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Proof. Coupling argument. □ 

Lemma 6.2. Let p > p c , and v be a probability measure on R + satisfiying the 
assumptions of the shape theorem. Denote by f m ; n the essential infimum of v, 
Vmean the mean of v , by n ViV the norm associated to the first passage-percolation 
model with S v = in the random environment given by the infinite cluster of 

Bernoulli(p) percolation, and by fi p the norm associated to the chemical distance 
for Bernoulli(p) percolation. Then: 

Proof. The lower bound is obtained by a coupling argument, the upper one by 
estimating the travel times on the minimal paths for the chemical distance. □ 

Theorem 6.3 (Flat edge result). Here, we work on 1? . Let p > p c = 1/2. Denote 
by p£ the critical threshold for oriented percolation on 1? . Suppose that S v = v® E , 
where v is a probability measure on R + such that f m ; n7 the essential infimum of 
v, is positive, satisfies the assumption of the shape theorem. Denote by A p the 
associated asymptotic shape, and by O(r) the ball with radius r for \\.\\i- 
1 



A p C O 



• If pv(vmiii) < Pc, then A p C int ( O 

• If pv{vmin) > pi, then 

A P n fr (o (— ) ) n (M+) 2 = 



l 



-M„ N q 



where 
q = pis(v min ), 

a q is the asymptotic speed in supercritical oriented percolation 
with pamameter q, 

" \2 ^2' 2 + ^2) ' 



denotes the segment line with extremities —^—M„, —^—N„. 

^min H ^min ^ 



Corollary 6.4 (Flat edge for the asymptotic shape for the chemical distance). 
Here, we work on Z 2 . Let p > p c = 1/2. Denote by pi the critical threshold for 
oriented percolation on 1? . Denote by A p the associated asymptotic shape, and by 
O the unit ball with radius 1 for \\.\\i. 



• A p C O. 



• If p < p c , then A p C int (O) . 

• If p>pl, then 

A p n/r(o)n(K+) 2 = [Mp,N p ], 

with the same notations as in the previous theorem. 
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The following pictures represent the set B n at time n — 450 for the chemical 
distance. Note the absence of flat edge for p = 0.55 < p c . 





p = 0.55 



p = 0.7 



Proof of theorem \6.'A Note that the first point is a direct consequence of lemma 
16.11 by comparing our model with the classical case and passage times with law 

Note then, by lemma 16.11 that the asymptotic shape of the first-passage perco- 
lation for some v®^ on the infinite cluster of Bernoulli(p) percolation is included 
in the asymptotic shape of classical first-passage percolation for the same v® m . 
Thus the second point is a consequence of theorem 1.3 (ii) in Mar02 , which is a 
consequence of the strict inequalities for the time constant in classical first-passage 
percolation established by van den Berg and Kesten in |vdB K93 . 

For the third point, let us examinate first the proof of the inclusion 



A p n fr 



D 



Yd 



-M. 



1 



-N n 



As noted in Kcs86 in the proof of theorem 6.13, the proof of the analogous result 
in the special case of Richardson's model treated in |DL81| lemma 6-13 by Durrett 
and Liggett can immediately be adapted. 

It remains now to prove the exact length of the flat edge by proving: 



A p n fr 



D 



n (R+) 2 c 



-M. 



But once again, the stochastic comparison lemma and the exact determination 
result in the classical case - theorem 1.3 in |Mar02j - give the announced result. □ 



7. Examples 

7.1. Exponential passage times. This subsection is devoted to the special case 
where the passage times are independent variables with an exponential distribution 
£(A). When p = 1, we recover a Richardson model, which is one of the first studied 
examples for first-passage percolation - see [Ric73 . 

By theorem 14. II /i is a norm for S u — (£ (A))® E and each p > p c . Let us denote 
by A(p : A) the asymptotic form corresponding to = (£(\))® Ed ■ By the scaling 
properties of the exponential laws, it is easy to see that A(p, A) = XA(p, 1) holds 
for each A > and p > p c . Let us prove that A(p, 1) C A(l,p) = pA(l, 1). 



CHEMICAL DISTANCE AND FIRST-PASSAGE PERCOLATION 



25 



The main tool is an appropriate coupling. Let (r] e ) e( z E d, i x e) e eE d i { z e) e £E d be 
independent random variables with for each e G E rf : r\ e ~ £(p), x e ~ £(1 — p) and 
z e ~ £(1). Now define 

= l{77e<a e } ? 7e + l{l) e >x e } Z e 

and 

w e = l{j7 e <a:e}- 

On one hand, 

P(T]' e >t,W e = l) = P(t] e >t,T) e < X e ) 

/+oo 
pe-™P(x e > y) d\(y) 

/+oo 
pe -py e -(i-p)y d\(y) 

= pe" 1 . 

Particularly, P(uj e = 1) = P(f]' e > t, ui e = 1) = p. On the other hand 

P(tf e >t,LJ e = 0) = P(z e >t,T) e < x e ) = P{z e > t)P{ Ve < x e ) = e _ *(l - t) 

Thus, rj' e and u> e are independent random variables, with rj e ~ £ (1) and w e ~ Ber(p). 

Consider now the set of wet points B t and B' t respectively associated to ((1, r] e ) ee ^d) 
and {{<jJ e , 7 le))eeE d - Since rf e — r\ e when the bond e is open, it follows that 

B'(t) C B(t) 

always holds. 

Dividing by t, and letting t tend to infinity, we get A(p, 1) C A(l,p) = pA(l, 1). 

7.2. Dependent \ 2 passage times. We then give an example of unbounded de- 
pendent passage times which lead to an asymptotic shape. It is also an example in 
which our theorem 14 . 31 makes us able to prove that /i is effectively a norm. 

Theorem 7.1. Let (^Q) ra ez d ,ie{i,....d} be a centered Gaussian process such that 

Vi,je{i,...,d} Vfc,neZ d El'4 = lljl| 

and 

£ Yl |EX^|<+oo. 

/cgZ d l<i<j<d 

To avoid degenerate cases, we also suppose that 

Vie{l,...,d} E pf*) 2 >0. 

./Vow consider the passage time (^ e )eeE d defined by r}ik,k+ei} = (X|,) 2 . J/ u;e denote 
by S v the law of (r/ e ) e ^d, then S v satisfies to the assumption of theorem 15., VI for 
each p > p c , so there is convergence to the asymptotic shape for the associated 
growth model. 

Proof. It is immediate that S v is shift-invariant. Since the covariance is summable, 
it tends to at infinity; it follows that S v is mixing, hence it is ergodic - this result 
is attributed to Maruyama and Fomin by Sinai' |Sin76j . 

Let us first recall some useful estimates whose proofs can be found in |Oarf)2| . 
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Lemma 7.2. Let X be a R n -valued centered Gaussian vector with covariance ma- 
trix C. Let us denote by p(C) the spectral radius i.e. the largest norm of an 
eigenvalue of C . Then, for each a 2 > p(C), we have: 

(7) P{\\X\\\ > no 2 ) < e - nh{ ^\ 

where h(x) = \{x — lnx — 1). The function h is increasing and positive on (1, +oo), 
with +oo as limit at +oo. 

Lemma 7.3. Let X be a n-dimensional Gaussian vector with positive definite co- 
variance matrix C . Let us denote by T(C) the spectral gap i.e. the smallest eigen- 
value of C . Then, for each a 2 < T(C), we have 

(8) P{\\X\\l < no 2 ) < e - nh( ^\ 

where h(x) = ^(x — lux — 1). The function h is positive and decreasing on (0, 1), 
with an infinite limit at 0. 

We will now prove that (?7e) e eE d satisfy to the assumption of theorem ED 
Let us for convenience identify E d with Z d x {1, . . . , d}: the point (x, i) 6 Z d x 
{1, . . . d} will be identified with {x, x + e{\ £ E d . We define 

S = E E |E XqX^ I < +oo and a 2 = inf {E (X^) 2 > 0; i £ {1, . . . , d}}. 

fceZ d l<i<j<d 

Let A C Z d x {1, . . . d} and e > 0. We must bound 

sJ v en s; J2ve<£\M)=P\ E iKf<^\ 

V eGA / \(k,i)eA 

Since E E |E X^Xft < ^ ^ |E *o^|, we can find N £ Z + 

keNZ d \{0} l<i<j<d k>N l<i<j<d 

such that 

E E |EJMl<y. 

k£NZ d \{0} l<i<j<d 

For each (k, i) £ {0, N- l} d X {1, . . . , d], we can define A kil = An ((k + NZ d ) x {i}). 

By the pigeon-hole principle, there exists (fcoj^o) such that |^4fc ,i | > jnJs- Let X 
be the |^4fc 0; i |-dimensional Gaussian vector composed by the (X e ) e£ A ko tQ ] it is 
obvious that 

E (^) 2 <|A| £ ) < F(\\X\\ 2 2 <\A\e) 

v (A,t)6A / 

< ¥(\\X\\ 2 2 <\A ko , io \dN d e). 

By lemma 1731 

dN d e < T(C) =► PdllH 2 . < |A fe0 , l0 |diV d e ) < exp (-\A koM \h (f^)) , 
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where the covariance matrix for X is C — (E X\X\ )((k,i),(l,j))eA ko iQ xA ko 4 ■ But 

Z i E x i x i° 



T(C) > inf E(X*°) 2 - £ |EX*X*°' 

°' (i.io)eA fco , lo ;(i,t )?4(fc,<o) 



> a 2 - sup £ 

(fe,Jo)eA feo , io (; !l0 ) e A fc0 , i0 ;(;, l0 )^(fc,») 

> a 2 - sup z iE^r-zi 

(fe, t0 )e J 4 feo , io ( Mo)6j 4 fco>lo ;(Mo)^(M 

> a 2 - £ |EXj>J#| 



> a 2 



k£NZ d \{0} 

2 ~ 2 ' 



It follows that 



(9) Ve e ( 0l ^) P £ W '<|A| E <e XP (-^(|)|A|). 

V 7 \(fc,i)£A / V 7 

We can now claim that the assumptions of theorem 14.31 are fulfilled. Since 
lim^^o h[x) = +oo, we have Kq — 0. Then, the application /i associated to p 
and S v is a norm for each p > p c . 

With the help of lemma VHH it is not difficult to see that 



eGA / 



S„ U e > n e > |A|a 2 < e 



-|A|M- 



holds for each a > S. It follows that (H a ) holds for each a > 1, and hence that the 
assumptions of theorem 15. 31 are fulfilled. □ 

7.3. Cost ans speed for road networks. Imagine that there exist n companies 
(Ci)i<i<n that want to build roads on the same Z d lattice: Company Ci, indepen- 
dentlty from the others, tries to build a road on edge e with probability p t (set then 
^ei = 1) an d leaves the edge e empty with probability 1 — p.; (set then to e i = 0), 
all constructions being independent. When two or more companies want to build a 
road on the same edge, they share the possession of the road. A client of company 
Ci can only use roads at least partially own by C,-. Suppose that each pi > p c (d) to 
ensure that almost surely, company d has a real (infinite) road network. Portions 
of roads own by comp d that are not linked to the infinite cluster cannot be used 
by its clients. Set, if e = {a, b} G E d is an edge: 

n 

Ve = E^e,il|Ci(o)| + |Ci(6)|=+oo- 
i=l 

Then ?y e counts the number of companies that offer the edge e as a possible road 
for their clients. Now, associate to the edge e a random value: 

Ve.i = f(i,Ve), 

where / is a deterministic function. 

For instance, we can imagine that f(i, r/ e ) represents the costs for company Ci to 
build a road on the edge e: this cost depend on the specificity of company Ci (that 
is why / depends on i), and also on the total number rj e of companies that want to 
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build a raod on th edge e (their can share the costs of construction for instance). 
In this case, the function / is decreasing in the second variable. 

Another example can be the following: f(i,r] e ) represents the time needed for 
clients of company d to cross the edge e: the larger the number of companies that 
own this edge is, and the larger the number of users to take this portion will be, 
and the longer the time needed to cross the edge will be. In this case, the function 
/ does not depend on i and is increasing in the second variable. 

The next result gives the existence of an asymptotic shape in this context. In 
other words, the travals times for a client of d will give raise to a deterministic 
asymptotic shape that does not on the random environment (neither his company 
road network, nor the other networks), but only on the densities (pi)i ie i<n of the 
distinct networks and on the function /. The same result is of course available for 
the cost point of view. 

Theorem 7.4. Let fl = (fls) n and P = P Pl x P P2 • • • x P Pn , with p { > p c for each 
i G {1, . . . , n} and f : {1, . . . , n} x {0, . . . , n} — ► M + with f(i, n) > as soon as 
n 0. We say that two points k and I of 7L d are i- connected if there is a path 
between k and I such that w e ,j = 1 for each bond e which is used in the path. For 
x G Z d , we denote by Ci{x) the i-connected component of x. Let us also define for 
e = {a, b) G E d : 



As before, for (ui, rf) G fl, and (x, y) G Z d xZ d , we define the travel time di(x, y)(uJ, rf) 
to be 



where the infimum is taken on the set of paths whose extremities are x and y and 
that are i-open in the configuration w. For t > and i G {1, . . . , d}, we note 



n 




i=l 



infV /(«,%), 



B\ = {ke 'L d ;d l {Q 1 k) < t}. 
For i G {1, . . . , d}, let us note by Pi the probability measure defined by 




t^+oo t 



where the convergence holds for Hausdorff 's topology. 
Proof. Fix i G {1, . . . , n}. For e = {a, b} G E d , let us define 




e£7 
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where the infimum is taken on the set of paths whose extremities are x and y and 
that are «-open in the configuration lu. For t > 0, we note 

B* t ={ke Z d ;d*{0,k) < t}. 

It is clear that 

Pi(Vf>0 B\=B* t ) = l. 

Since the passage time {vit)esM d are independent from the environment {^k,i)kez d i 
it will be possible to use the preceding theorems. 

The first step is to check that the law of (rj*) eeK d is shift-invariant and ergodic 
under the translations of 1 d . Obviously 

® P w = ( ® (1 ~ P 3 )5o + PA)®®* 

r II n) \>\ ,J V je{l,...,n}\{i] 1 ° ' 

is shift-invariant and ergodic under the translations of % d . For simplicity, we note 
L=( <g> 8 +p j 6 1 ) mi . We note 

j'6{l,...,n}\{t} 

F(UJ) = (f(i,l+ W e jl|C J (a( e ))| + |C 3 (b(e))l=+oo)) eeE d, 

i6{l,...,n}\{i} 

with the notation e = {a(e), 6(e)}. 

The law of {j)e) e eM d is precisely the image of L by the transformation F. Let us 
denote by this law. It is easy to sec that F o 9 U = 6 U o F holds for each u 6 Z d , 
so that we have the commuting diagram: 



I' I' 

(R+f — ^ (K+f 
For any event A £ S((R+) Ed ), we have 

= L^oF)- 1 ^)) 
= Lfffo^)- 1 ^)) 

= l^-Hf- 1 ^)) 

= L(F- 1 (A)) = 5,(A). 

Then, is invariant under the translations. 

It is now clear that the dynamical system ((R + ) E , S((K + ) E ), S v , 9 U ) is a factor 
of the dynamical system (({0, l}™" 1 )^, S(({0, 1}™- 1 ) 1 "), L, 6 U ). Since the dynam- 
ical system (({0, l}"" 1 )^, S(({0, l}"" 1 )^), L, 6 U ) is ergodic, it follows that the 
dynamical system ((K+) E<i , i3((R+) E<i ), S v , U ) is ergodic too. 

Now, since min/(i, 1 + .) > and max f(i, 1 + .) < +oo, the conclusion follows 
from the remark that we have done before the proof of theorem 15.31 □ 
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